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Abstract 

A fourth generation of Standard Model (SM) fermions is usually considered unlikely due to 
constraints from direct searches, electroweak precision measurements, and perturbative unitarity. 
We show that fermion mass textures consistent with all constraints may be obtained naturally 
in a model with four generations constructed from intersecting D6 branes on a T^/(Z2 x Z2) 
orientifold. The Yukawa matrices of the model are rank 2, so that only the third- and fourth- 
generation fermions obtain masses at the trilinear level. The first two generations obtain masses 
via higher-order couplings and are therefore naturally lighter. In addition, we find that the third 
and fourth generation automatically split in mass, but do not mix at leading order. Furthermore, 
the SM gauge couplings automatically unify at the string scale, and all the hidden-sector gauge 
groups become confining in the range 10^^-10^^ GeV, so that the model becomes effectively a 
four-generation MSSM at low energies. 
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I. INTRODUCTION 

Current high-energy experimental data supports just three generations of chiral fermions. How- 
ever, the existence of a fourth generation remains viable provided the mass of the fourth-generation 
neutrino v' satisfies rriyi > \Mz-, and the corresponding charged fermion masses m^', m^', and m,-/ 
lie in the correct mass ranges to avoid constraints from direct searches and precision electroweak 
measurements Experimentally, various anomalies indicate possible manifestations of a fourth 
generation. Recent analysis from WMAP7 |2] suggests a number of relativistic neutrino flavors 
that is greater than three, A'efr = 4.34j'gg. Explanations have been proposed based on sterile 
neutrinos with masses at sub-eV scales [3| (see also Refs. for earlier discussion of astrophysical 
implications of the fourth generation). In addition, recent observations from MiniBooNE 5| and 
the observation of a reactor antineutrino anomaly may indicate the existence of one or more 
light, sterile neutrinos. Furthermore, the existence of a fourth generation generically introduces 
additional sources of CP-violating effects U|i^ which are important for baryogenesis and could 
explain observed anomalies in flavor physics [10,]. 

Although the existence of a fourth generation is still viable, it has long been considered unlikely. 
In part, this predisposition relies upon the perceived unnaturalness of the allowed fourth-generation 
fermion masses, which are strongly constrained by the limits provided by direct searches and pre- 
cision electroweak measurements. In particular, these constraints mandate a large v' mass in com- 
parison to the tiny masses of the observed three neutrinos. Moreover, the fourth-generation quark 



masses are constrained by the perturbative unitarity of heavy-fermion scattering amplitudes [111] to 
be < 500-600 GeV. Thus, the t' quark mass m^' can only lie between ~ 2-3 times the m^, a small 
difference when compared to the mass hierarchy m(/mc~ 135. Indeed, for the known families of 
quarks and leptons, progressively larger mass splittings between each sequential family appears to 
be the rule, a pattern that cannot hold in four-generation models due to the limited allowed range 
for rriti . Furthermore, contributions to the oblique S and T parameters from an additional chiral 
fermion generation are intolerably large unless the masses of the t' and h' quarks are equal to within 
rat' — ray ~ (1 -|- ^Inyj^) x 50 GeV, where ra^ is the Higgs mass [l^. Moreover, mixing between 
the third- and fourth-generation quarks is often presumed small to accommodate the unitarity of 
the CKM matrix. These features of four-family models appear to require significant fine tuning to 
explain, and therefore a fourth generation of chiral fermions is usually deemed improbable. 

In fact, the pioneering work of Ref. [3] showed that fourth-generation fermions are phenomeno- 
logically perfectly viable, apart from mild tuning of the masses within the fermion 8X1(2)^;, doublets. 
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The direct search hmits quoted therein are 111^1,111^1 > 100 GeV and niti^my > 258 GeV, and the 
electroweak precision-based constraints are given by m^i — / ~ 30-60 GeV and the m^/ — my 
sphtting mentioned above, while the most stringent quark-mixing bounds are actuahy those with 



the first generation, |V^fr 
tailed calculations {e.g., 
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{iLd\ ^ 0.04. Subsequent bounds from both observation and more de- 



141]) modify these values slightly, but the possibility of a fourth chiral 



fermion generation remains vibrant and awaits the decisive verdict of the Large Hadron Collider. 

Currently, string theory is the most promising approach for a unification of gravity with quantum 
mechanics. As such, it should provide a first-principles explanation of the detailed properties of our 
universe. However, as is well known, string theory exhibits a huge degeneracy of vacua, resulting 
in the so-called string landscape. For many string compactifications, replication of chiral matter 
seems to be a generic feature. At present, the principle that sets the precise number of observed 
chiral generations remains unknown. Indeed, many consistent string vacua can be constructed 
that closely resemble the Minimal Supersymmetric Standard Model (MSSM) but contain different 
numbers of generations [15,]. 



One particu 
reviews, see 



arly promising class of string compactifications employs D branes on orientifolds (for 



16l4l9l|). In such models chiral fermions — an intrinsic feature of the Standard Model 



(SMV— arise from configurations with D branes located at transversal or bifqld / conifold singulari^ 



ties 



20| and strings stretching between D branes intersecting at angles 
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picture, with magnetized D branes 



22l | (or, in its T-dual 



23l425l| ) . Within the framework of D-brane model building on 



toroidal orientifolds, the worldsheet areas Aabc spanning the region bounded by D branes (labeled 
by a, 6, c) that support fermions and Higgses at their intersections give rise to Yukawa couplings 



Yabc ~ exp(-^abc) 



22 



2a]. This pattern naturally encodes the hierarchy of Yukawa couplings. 



Indeed, models exist in the literature in which realistic mass matrices and CKM mixings for both 
three- and four- generation models may be obtained [3]. While accommodating full rank-3 or rank- 
4 Yukawa matrices is possible, one may nevertheless question whether the degree of fine-tuning 
of parameters required is too large to provide a completely natural explanation of the observed 
fermion mass and mixing hierarchy. 

In this paper we construct a four-generation MSSM with intersecting D6 branes on a 
r^/(Z2 X Z2) orientifold. We show that the Yukawa matrices are rank 2, and thus only the 
heavier two generations receive masses at the trilinear level, while the lighter two generations can 
receive masses from higher-order couplings.^ Moreover, the third- and fourth-generation fermions 



^ In some cases, use of the term "naturally lighter" must be carefully considered: Note, for example, the second- to 
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automatically split in mass, but we find that they do not mix at the trilinear level and therefore 
do not generate 0(1) mixing angles. Thus, one can obtain fermion mass textures naturally consis- 
tent with fourth-generation constraints. Finally, the tree-level gauge couplings automatically unify 
at the string scale, and the hidden-sector gauge groups introduced to satisfy tadpole constraints 
become confining in the range 

10i3_io^6 QgY^ leaving only the MSSM at low energies. 



II. A FOUR-FAMILY MSSM 
A. Model Building on T^/{'L2 x Z2) 

Let us begin by considering Type IIA string theory compactified on a T^/(Z2 x Z2) orien- 
281 ]. where the six-torus is factorizable in terms of two-tori as = X X T^, and 



tifold [27 

the i^^ two-torus, i = 1,2,3, has complex coordinates Zj. The 9 and w generators for the orbifold 
group Z2 X Z2 act on the complex coordinates of as 

■■ {Z1,Z2,Z3) ^ {-Zl,-Z2,Z3) , 

CO : {zi,Z2,Z3) ^ {zi,-Z2,-Z3) . (1) 

The orientifold projection is implemented by gauging the symmetry QR, where is worldsheet 
parity, and R is given by conjugation: 

R : {zi,Z2,Z3) ^ (zi,Z2,Z3) . (2) 

As a result, one finds four kinds of orientifold 6-planes (06 planes), corresponding to the actions 
QR, 0,R9, QRlo, and VlROuj. Only two kinds of compL 



projection for a two-torus: rectangular and tilted 



ex structures are consistent with orientifold 



28|, with A = 0, 1 for an untilted or tilted i^^ 



torus, respectively. In the following, we only consider rectangular two-tori (but retain parameters 
allowing for tilt to permit the generalization of the model presented here). The homology three- 
cycles for a stack a of D6 branes and its orientifold image a' are given by 

m = n«[«^] + ^'^%M, pa'] = nK[«^] - ^~^%M . (3) 

i=l i=l 

The Ramond-Ramond (R-R) tadpole cancellation requires the total homology-cycle 
charge of D6 branes and 06 planes to vanish [2^ , which may be as expressed as 

J2 NaiHa] + J2 ^'^Pa'] " 4^06] = . (4) 



third-generation mass ratio mc/rriT ~ 0.7. However, a strong hierarchy remains when one compares only fermions 
of the same type, e.g., m^<^ mr- 
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It is convenient to use the parameters 



Aa 



''a'-a''a i 



1/2,3 



^a. ~ ''a"'a''a ) ~ '•a''a^a ) 



Da — '^a^a^a i 
1,2,3 



(5) 



which were first introduced in 



291 ] . The tadpole conditions can then be expressed as 



a a a a 

Completely eliminating the R-R tadpoles actually requires the cancellation of D-brane 



as classified by K theory, which for tori of arbitrary tilt for the Z2 x Z2 orientifold requires 
J^iVaiaG4Z, ^Ar„5ae4Z, ^iVaae4Z, ^iVa^ae4Z. 



13, 



(6) 



charges 



30, 



3li 



(7) 



In order to preserve J\f = 1 supersymmetry in four dimensions, the rotation angle of any D brane 
with respect to the orientifold plane must be an element of SU(3) [2l|, |28(. Consider the angles 6^ 
between each brane a and the i?-invariant axis of the i^^ torus; one requires that 9^^ + 6'^ + 6^ = 
mod 27r, which means sin(6'^ + + ^a) = and cos(6'^ + ^'a + ^a) > 0- Define 



tan ( 



(8) 



where R\ 2 are the radii of the i torus and Xi '- 
supersymmetry conditions can then be recast as 



are the complex-structure moduli. The 



29| 



XAAa + XsBa + XcCa + XoDa = , 
Aa/XA + Ba/XB + Ca/xc + Da/xD < ; 



(9) 



where xa, xb, xc, xd are complex-structure parameters, which are given by 

XA = X, = A • 2^2+/33/^2X3, xc = A • 2^i+''VxiX3, xd = X- 2^^+^^ /xiX2 , (10) 

where A is a positive parameter introduced [29] to put all the variables A,B,C,D on an equal 
footing. However, among the Xj, only three are independent. 

The initial U(A'^ri) gauge group supported by a stack of Na identical D6 branes is broken down by 
the Z2 X Z2 symmetry to a subgroup U(iVa/2) [28] . Chiral matter fields are formed from open strings 



whose two ends attach to different stacks. By using the Grassmann algebra [oj] [bj 



-[bj][ai 



and [oj][aj] = — [&j][^j] = 0, one can calculate the intersection numbers between stacks a and b and 
obtain the multiplicity (M) of the corresponding bifundamental representation: 

3 

= la, = [Ua] m = 2-'^ Hiny, - niQ , (11) 



M 



Na Nb 



2 ' 2 



i=l 
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where fc = /3i+/?2+/33 is the total number of tilted tori. Likewise, stack a paired with the orientifold 
image b' of b yields 

M ^) = lay = [n„][n,,] = -2-^ + tATo) ■ (12) 

^ ^ i=i 

Strings between a brane in stack a and its orientifold image a' yield chiral matter in the anti- 
symmetric and symmetric 05 representations of the group \]{Na/'2,), with multiplicities 



M[(aA,a)L\ = ^Ia06, M[{aA,a +as,a)L] = ^ (^ha' " ^hosj , 

SO that the net numbers of antisymmetric and symmetric representations are given by 

M{aA,a) = I (laa' + llaOb) = "2^"'= [(2 A " l)ia -Ba-Ca- Da] , 



(13) 



2 V 2 

M{as,a) = \ (laa' " ^^aOe) = -2'-''[{2Aa + l)Aa + Ba + Ca + Da] , (14) 

where 

3 

Iaa' = [na][Ua,]=-2^-'^llnili, (15) 

1=1 

Ia06 = [Ha] [Hoe] = 2'-'=(ia + Ba + Ca + Da) ■ (16) 



Note that the expressions Eqs. (jl3p and (jl4p indicate the presence of non-chiral pairs of matter in 
the antisymmetric representation that are masked by Eq. ()14p . which gives only the net number. 

A zero intersection number between two branes implies that the branes are parallel on at least 
one torus. At such types of intersection, additional non-chiral (vectorlike) multiplet pairs from 
ab + ba, ab' + b'a, and aa' + a' a can arise |3^^. The multiplicity of these non-chiral multiplet pairs 
is given by the remainder of the intersection product, neglecting the null sector. For example, if 
Kll - nlll) = in la, = [nj [n^] = 2-fc YiUiKll - then 

3 



M 



Na N,\ (Na iVAl fr, 

1~'T/ l^'T/ = V\y^<^''b - ^bL) ■ (17) 

^ ^ / -I i=2 

This result is useful since one can fill the spectrum with this matter without affecting the required 
global conditions, because the total effect of the pairs is zero. Typically in this type of model, the 
Higgs fields arises from this non-chiral matter. However, as we shall see for the model constructed 
here, the Higgs fields appear in the chiral sector of the model. 



^ Representations (Antia + Antia) occur at intersection of a with a' if they are paraUel on at least one torus. 
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The total non-Abelian anomaly in intersecting brane-world models cancels automatically when 
the R-R tadpole conditions are satisfied, but additional mixed anomalies may be present. For 
example, the mixed gravitational anomalies that generate massive fields are not trivially zero [28|]. 
These anomalies are cancelled by a generalized Green-Schwarz (G-S) mechanism that involves 
untwisted R-R forms. The couplings of the four untwisted R-R forms B2 to the U(l) field strength 
Fa of each stack a are given by 

NJlnlnl [ Bl A trF^, Naulllnl [ S| A trF„ , 
Jam Jm4 

Naulnlll I Bl A trF,, -NJlllll [ BI A trF, . (18) 

These couplings then determine the exact linear combinations of U(l) gauge bosons that acquire 
string-scale masses via the G-S mechanism. If U(l)x is a linear combination of the U(l)s from 
each stack, 

= J^CaU(l)a, (19) 

a 

then the corresponding field strength must be orthogonal to those that acquire G-S mass. Thus, if 
a linear combination U(l)x satisfies 

Y,CaNaBa = Q, Y,CaNaCa = 0, 

a a 

Y,CaNaDa = 0, ^ C„iVa^ = , (20) 
a a 

the gauge boson of U(l)x acquires no G-S mass and is anomaly-free. 



B. The Model 



We present the D6-brane configurations, intersection numbers, and complex-structure parame- 
ters of the model in Table HI and the resulting spectrum in Table HH With this configuration, all 
R-R tadpoles are cancelled, M = 1 supersymmetry is preserved, and K-theory conditions are satis- 
fied. The resulting model has gauge symmetry [U(4)c x U(2)l x U (2) r] observable x [USp(8)3]hidden, 
a left-right symmetric Pati-Salam symmetry that provides (as seen below) a convenient origin for 
the SM hypercharge U(l)y. 

In this model the anomalies from the three global U(l)s of U(4)c, \J{2)l, and U(2)i^ are each 
cancelled by the G-S mechanism, and the gauge fields of these U(l)s obtain masses via the linear 
B A F couplings, except for one linear combination U(l)x = U(l)c + 2[U(1)2l + U(1)2_r]. The 
effective gauge symmetry of the observable sector is then SU(4)c x S\J{2)l x SU(2)/j x U(1)x- 
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The Pati-Salam gauge symmetry must then be broken to the SM. The breaking is accomphshed 
by brane sphtting, which corresponds in the field theory description to giving VEVs along D-flat 
and F-flat directions to some of the adjoint scalars that arise via open-string moduli associated 
with each stack of D branes |33l ]. The a stack of D6 branes may be split on any of the three two-tori 
into al (quark) and a2 (lepton) stacks with Nai = 6 and Na2 = 2 D6 branes, respectively. In same 
manner, the c stack of D6 branes is split into stacks cl {I^l = +| fermions) and c2 {I^l = 
fermions) such that Nd = 2 and Nc2 = 2, respectively. After the splitting, the observable gauge 
symmetry is broken to SU(3)c x SU(2)j^ x U(l)/3^ x U(1)b-l x U(l)x) where the gauge bosons of 
U(l)/3^, \J{1)b-l, and U(l)x remain massless. The U(l)/3^ x \J(1) b-l x U(1)x gauge symmetry 
may then be broken spontaneously if, for example, vectorlike particles with the quantum numbers 
(1, 1, i, -1, 3) and (1, 1, -i, 1, -3) under the SU(3)c x SU(2)l x U(l)/3^ x \J{1)b-l x U(1)x gauge 
symmetry from the 0202 intersections $ in Table HI]) receive VEVs,^ leaving only the single 
U(l) corresponding to SM hypercharge, 

U(l)y = ^[U(l)al - 3U(l)a2 - 3U(l)cl + 3U(l)e2] , (21) 

and the overah gauge group SU(3)c x SU(2)l x U(l)y x USp(8)3. 

The observable sector of the model then possesses the gauge symmetry and matter content of 
a four-generation MSSM with an extended Higgs sector. Note that the four pairs of Higgs fields 
in this model appear in the chiral spectrum, rather than arising as vectorlike matter from N = 2 
subsectors, as is the case for the models studied in [34], as well as most models in the literature. 
The extra matter in this model includes fields charged under the hidden-sector gauge groups, 
vectorlike matter between nonintersecting pairs of branes, and the chiral adjoints associated with 
each stack of branes. In addition, matter occurs in the symmetric triplet and antisymmetric singlet 
representations of 811(2)^. 

In order for the model to reproduce just the MSSM at low energies, the gauge couplings must 
unify, and all extra matter besides the MSSM states must become massive at high-energy scales. 
Furthermore, since the MSSM contains just one pair Hu^d of light Higgs doublets, the Higgs po- 
tential must be fine tuned. Since the four Higgs fields in this model appear in the chiral sector 
rather than as = 2 vector like multiplets, no intrinsic /i parameter (whose real part corresponds 
geometrically to the separation between stacks h and c) occurs since the b and c stacks intersect 
on all three two-tori. A single pair of light Higgs doublets may be obtained by fine-tuning the ^ 



^ The a2c'i intersections were incorrectly indicated for this purpose in Ref. [l^. 
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TABLE I: D6-brane configurations and intersection numbers for a four-famiiy Pati-Saiam modei on a Type- 
IIA r^/(Z2 X Z2) orientifoid, witlr no tilted tori. Tlie complete gauge symmetry is [U(4)c x U(2)i x 
U(2);j]obscrvabic X [USp(8)'^]hiddcn, and the complex-structure parameters are xi — 2x2 = 2x3- The parame- 
ters Pi are the /3-function coefficients for the USp(8)j gauge groups. 







U(4)c X \J{2)l X U(2)« X 


USp(8)^ 












N 


(n^,m^) X {ii?,'w?) X (n^,m^) 


ns 


Ua 


6 


6' 


c 


c' 


1 


2 


3 


a 


8 


(0,-1) X (1,1) X (1,1) 








4 





-4 





1 


-1 





b 


4 


(2,1) X (1,-1) X (1,0) 


2 


-2 






4 








1 


-2 


c 


4 


(2,-1) X (0,1) X (1,-1) 


-2 


2 










-1 





2 


1 


8 


(1,0) X (1,0) X (1,0) 




Xi = 


= 2, X2 = 


1 


X3 


= 1 




2 


8 


(1,0) X (0,-1) X (0,1) 




/3 


9 _ 

1 - 




-3, 






3 




3 


8 


(0,-1) X (1,0) X (0,1) 








P 


9 _ 


-2 







term, which may be generated in the superpotential W via the higher-dimensional operators 

ijkl 

W D y^^SlSi^HM , (22) 

where y]^^'' are Yukawa couplings, and M^t is the string scale. In this scenario, the singlets S'^ 
may obtain string or GUT-scale VEVs while preserving the D-flatness of U(l)/j, and the singlets 
S*^ may obtain TeV-scale VEVs while preserving the D-flatness of U(l)i^ (and hence are tied to 
electroweak symmetry breaking). The precise linear combinations that produce the two light Higgs 
eigenstates Hu,d obtained by fine-tuning the Higgs potential via Eq. ([22]) are then correlated with 
the pattern of Higgs VEVs necessary to obtain Yukawa matrices for the quarks and leptons, 

where v^^^ = (^Hl^y 

The gauge coupling constant associated with a stack P is given by 

9m, = |Re(/p)|, (24) 

where fp is the holomorphic kinetic gauge function associated with stack P, given by 13] 

fp = (np UpUpS - Up Ip Ip - Up Ip Ip - Up Ip Ip u^) , (25) 

with Kp = 1 for SU(A''p) and Kp = 2 for USp(2A''p) or S0(2A''p) gauge groups, and s and u are 
moduli in the supergravity basis. The holomorphic gauge kinetic function associated with the SM 
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TABLE II: The chiral and vectorlike superfields, their multipHcities and quantum numbers under the gauge 
symmetry SU(4)c x SU(2)l x SU(2)fl, x USp(8)i x USp(8)2 x USp(8)3. 





Mult. 


Qujiiituiri NuuibBr 


Oa 




^ ZH 


Field 


ab 


4 


(4, 2, 1, 1, 1, 1) 


1 


— 1 





Fl{Ql, Ll) 


ac 


4 


(4, 1, 2, 1, 1, 1) 


— 1 





1 


Fr{Qr, Lr) 


be 


4 


(1, 2, 2, 1, 1, 1) 





1 


— 1 


HI, 


al 


1 


(4, 1, 1, 8, 1, 1) 


1 








Xal 


a2 


1 


(4, 1, 1, 1, 8, 1) 


— 1 








Xa2 


ol 


i 


(i, z, i, i, 0, i j 


o 

u 


1 

i 


U 




63 


2 


(1,2, 1, 1, 1, 8) 





_x 





xi 


cl 


I 


(1, 1,2, 8, 1, 1) 








_1 


X 1 


c3 


2 


(1, 1, 2, 1, 1, 8) 








1 






2 


(1, 3, 1, 1, 1, 1) 





2 





Tj 

L 


bA 


2 


(1,T, 1, 1, 1, 1) 





-2 





Si 




2 


(1, 1, 3, 1, 1, 1) 








-2 


Ti 

^ R 


CA 


2 


(1, 1, 1, 1, 1, 1) 








2 


Sh 

^ H 


ao 


o 
2 


(4, z, i, i, i, i j 


1 

i 


1 

i 


U 






o 
z 


(^4, z, i, i, i, 1} 


— i 


1 

— i 


U 




1 

ac 


Z 


(4, i, Z, i, i, L) 


1 

i 


u 


i 


<Pi 




2 


(4, 1,2, 1, 1, 1) 


-1 





-1 




bb' 


2 


(1,T, 1, 1, 1, 1) 





-2 









2 


(1, 1, 1, 1, 1, 1) 





2 





Si 


cc' 


2 


(1, 1, 1, 1, 1, 1) 








2 






2 


(1, 1,T, 1, 1, 1) 








-2 




be' 


1 


(1,2,2,1,1,1) 





1 


1 


H' 




1 


(1,2,2, 1, 1, 1) 





-1 


-1 


h' 



hypercharge U(l)y is then given for this model by the combination 36|] 



fr = -/al + -fa2 + -/cl + -/c2 • 
D Z 2 Z 



(26) 



The complex-structure moduli [/* are obtained from the conditions to maintain A/" = 1 supersym- 
metry. Consistent with the complex-structure parameters xi = 2, X2,3 = l! this model has 



[/I = 2i , C/2 = i , C/3 



^ . 



(27) 



Note that the conditions for preserving = 1 super symmetry do not actually fix one complex 
structure modulus, , and we have chosen a value for it which results in gauge coupling unification. 
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In a completely realistic compactification, it would be necessary to find some mechanism to stabilize 
this modulus, as well as the Kahler and open-string moduli, and the dilaton. In the supergravity 
basis. 



Re(s) 
Re (n^ ) 




Vim C/i Im C/2 Im C/3 



ImW 



Im Ul" Im 



jjk jjl 



{j,k,l) = (1,2,3), 



(28) 



where (^4 is the four-dimensional dilaton, one finds 

1 e-<^4 I £-04 



Re(s) 



Re('u^) 



Re(n^ 



1 2e- 



Re('u3) 



1 26-"^^ 
2^ V2 



(29) 



27r V2 ' 2^ V2 ' 2^ V2 

Inserting these expressions into Eqs. (iM|) - (f25] ). one finds that the gauge couplings are unified 
at the string scale. 



5s =9l = lol = 2\/27re'^ 



(30) 



Choosing (/)4 = -2.295 so that e't'^ = 0.101, one obtains g^^^jj 0.895, the value of the unified gauge 
coupling at GUT scale for a four-generation MSSM. One further finds the string scale to be 



Mst = 7ri/2g<^4^pj ^ 4 35 y. iQiY Qgv, 



17 



(31) 



where Mpi is the reduced Planck scale, 2.44 x 10^^ GeV. 

Fixing the value of 4)4, also fixes the gauge couplings of the hidden-sector USp(8) groups at Mst: 



5? = (71 = 7.160, ff| = 3.580. 



(32) 



Using the /3-function coefficients shown in Table HI the renormalization-group equations (RGEs) 
for the USp(8)i^2 gauge groups produce strong coupling at a scale Ai^2 ~ 4.5 x 10^^ GeV, while 
the group USp(8)3 reaches strong coupling at A3 w 8.3 x 10^^ GeV. Thus, matter charged under 
these groups becomes confined into composite states and decouples not far below the GUT scale. 

Thus, at low energies this model is a four-generation MSSM, where the gauge couplings unify 
« 2.2 X 10^^ GeV. In addition, one finds matter charged under the hidden-sector gauge groups that 
becomes confined into massive bound states at high-energy scales. These matter representations 
all have fractional electric charges, Qem = particular, the fields charged under USp(8)3 

all have half-integer electric charges and confine into bound states with integer electromagnetic 
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40(] . If the lightest of these 



charge at a scale 0(10^^) GeV, similar to the "cryptons" studied in 
bound states is electromagnetically neutral and is metastable with a sufficiently long lifetime, it 
provides a good candidate for superheavy cold dark matter (CDM). This observation is particularly 
true for states with masses of order A3, since superheavy particles with masses in the range M = 
O(10^^~^'^) GeV might very well have been produced naturally through the interaction of the 



vacuum with the gravitational field during the post-inflationary reheating phase |41l] . 



III. YUKAWA COUPLINGS 



Formalism and Construction 



A complete form for the Yukawa couplings arising from D6 branes wrapping on a full compact 
space X X can be expressed as 



26 



Yijk oc i9 



r=l 



sir) 



where the proportionality constant cancels in mass ratios, and 

sir) 



1} 



Ar)) 



with r = 1, 2, 3 denoting the three two-tori. The input parameters are given by 



fir) 



rir) 



+ 



I. 



(r) 



+ 



be 



fir) Ar) Ar) 
^ab ^bc -"CI 



+ 



Tir)n{r) , Ar) Ar) Ar) Jr) 
dir) 



jM I Ar) Ar) Ar), 
.(r) _ l-'ab -'fee ^ca \ 



a' 



(dir)] 



(33) 



(34) 



(35) 



where the indices and k'^^^ label the intersections on the r^^ torus, = gcd[lj^J , 

Im), and the integer s^*"-* is a function of ji^\ and k^^^ corresponding to different ways of 
counting triplets of intersections. J^^-* is the Kahler modulus on the r^^ torus and a' is the string 

(r) 

, and ec correspond to the relative positions of stacks a, 6, 



(r) (r) 



' ^b 



tension. The shift parameters ea 
and c, while the parameters (j)a \ 4'b'\ and (p'^'^ are Wilson lines associated with these stacks. For 
simplicity and to exhibit the flexibility of our model, we set the Wilson lines to zero. The brane 
shifts and Wilson lines together comprise the open-string moduli, which must be stabilized in a 
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TABLE III: Index assignments labeling the intersections at which each chiral field is localized. 



Left-handed 


i = 


Right-handed 


j = 


Higgs field 


k = 


field 


{i(l),i(2),i(3)| 


field 








FliQlLl) 


{0,0,0} 


Fr{Qr, ifl) 


{0,0,0} 




{0,0,0} 


FliQhLl) 


{0,1,0} 


Fr{Qr, Lr) 


{0,0,1} 


H\Hl,Hl) 


{1,0,0} 


Fl{Ql,Ll) 


{1,0,0} 


fMqLl'r) 


{1,0,0} 


H\HlHl) 


{2,0,0} 


Ft{Qi,Li) 


{1,1,0} 


Fr^Qr, L%) 


{1,0,1} 


H\HtHi) 


{3,0,0} 



complete model by some mechanism. In this paper we treat them as free parameters; however, we 
discuss the stabilization of these moduli in the next subsection. For simplicity, let us define 

Mr) ( A"-) A"-) I rWpW I fir) (r). 
g(r) _ " y^ab + ^ca % + } 
Jr) j{r) j{r) " 
^ab ^bc -"ca 

Although these Yukawa coupling formulas originally refer to = T"^ x T'^ x T'^ , they may be 
generalized to /{'L2 x Z2) simply by including all of the orbifold images. However, in the present 
case the cycles wrapped by the orbifold images of a stack of D branes a are homologically identical 
to the original cycle wrapped by the stack a. In addition, the intersection numbers between the 
cycles defined on the orbifold turn out to be the same as the intersection numbers between those 
on the ambient torus. Thus, the above formulas for = x x may be used without 
change on T^/{Z2 x Z2). 

Since in this model the intersections are distributed over all three two-tori, they must be 
carefully labeled. Let us introduce the indices i = {A^\A'^\A^^^, j = J*-^-*}, and 

k = ^k^^\k^'^\k^^^^ to label the intersections at which the left-handed fields, right-handed fields, 
and Higgs fields, respectively, are localized. For this model, these indices can assume the values 

i(i)={0,l}, i(2)={o,l}, i{3)={0}, 

j« = {0,l}, J^'^ = {0}, j(3) = {0,l}, (37) 
fcW = {0,1,2,3}, A;(2) = {0}, k^^^ = {0} , 

as obtained by considering the intersection numbers /a?, ij,^^ on each torus. Then each 
left-handed, right-handed, and Higgs field is identified by a triplet of indices, as summarized in 
Table Hill The independent Yukawa couplings on each torus are then labeled: 

^ _ 1 . v(^) v(^) ^(1) ^(1) ^(1) Ci^) 

' ~ ^ • -'000' -'oil' -'^101' -'^110' -"^002' -'^013' -'^103' -'^112' V'-'°J 
r = 2: yi^Y^l ^2 ^ ^ol (39) 

r = 3: zi^Y^^l, Z2 = Fqi'd' (40) 
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being the only remaining components after applying the selection rules i^^^ +j^^^+k^^'^ = mod d^^^ , 



^(1) + + fc(i) = mod 2 , 
•(2) _^ j{2) _^ ^(2) ^ Q ^ ^ 

j(3) _^ ^-(3) _^ ^(3) ^ Q ^^Q^ ^ 



(41) 



The full Yukawa couplings are then given by a product of the couplings on each torus, 

3 



n 



r=l 



(42) 



Thus, the Yukawa matrices (including VEVs) in this model are of the form 



ijk^k 



YZ 
YZ 




(43) 



where YZ is the singular 2x2 matrix 



YZ = 



yizi yiZ2 
y2Zi y2Z2 



(44) 



and the values = are the VEVs of the Higgs fields. 

Since the same singular submatrix YZ appears in all sub-blocks of y, the full Yukawa matrices 
y are only rank 2, so that at most two eigenvalues of y can be different from zero. However, this 
result can be quite desirable, since the third- and fourth-generation quarks and leptons tend to be 
significantly heavier than those of the first and second generations. In principle, lifting the rank-2 
degeneracies (to induce small masses and mixings) can be achieved by using higher-order couplings 



than included here, such as the fourth-order terms considered in 
corrections 



431, 



4J], or supersymmetry loop 



45(1 . For example, the couplings 



„/ pi pk TTl 

yijk^L^L-^R-"- ■ 



(45) 



can perturb the Yukawa matrices to break the rank degeneracy, or additional trilinear couplings 
with the vectorlike Higgs fields H and H' can be induced by D-brane instantons. 

For this model, the 5 parameters that enter the -d functions of Eq. (j34p are given by 



^(1) 

^(2) 
^(3) 



2 



;(2) 



feW __ 
4 4 

b ^*^c J ) 



,(1) 



+ 



;{3) 



2^ 



(46) 
(47) 
(48) 
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For simplicity, as mentioned above we have taken c/)^^^ = (f)^'^^ = (f)^^^ = 0. Allowed choices for 
s^^^ are restricted by the requirement that the parameters S^^\i,j,k) must differ by an integer for 
triplets (i, j, k) forming triangles with the same area. In the present case, one finds two independent 
choices satisfying this constraint: s^^^ = i^^'^ and s^^^ = For definiteness, let us make the choice 
= jW. We then find 

"^000 ~ -^112 = ^0 5 ^on^ = ^101 = ^1 ' (^9) 

y(l) _ y(l) = y(l) _ y(l) = 

The Yukawa matrix then takes the simple form 

, (xoVo + X2V2)YZ {xiVi + X3V3)YZ 

y = \ I > (50) 

{xiVi + XsVs)YZ {X2V0 + XqV2)YZ 

where we again see that the Yukawa matrices are rank 2, so that two eigenvalues are zero. The 
four eigenvalues of y can be given in closed form as Ai = A2 = 0, and the two nonzero eigenvalues 

■^3 = ^(yi^i + y2^2)[(a;o + 2^2) (^^0 - V2) - A] , (51) 
-^4 = + y2^;2)[(a;o + X2){vq + V2) + A] , 

where 



A = V(xo - 3:2)^(^0 - V2Y + ^{xiVi + X3W3)2 , (52) 

from which one sees that a mass splitting between the third and fourth generations automatically 
occurs for generic nonzero values of the Higgs VEVs. In order to work with simple closed-form 
solutions, note that setting vq = V2 simplifies Eq. ([5T]) to just 

-^3 = ^(yi^i + ^2^2) [{xo + X2)vo - {xivi + xs^s)] , (53) 

A4 = ^(2/1^1 + ^2^2) [{xq + X2)V0 + {xiVl + X3V3)] . 

The ratio of the third- to fourth-generation fermion masses is then given by 

fu,d 



^ fu,d 



fu.d U,d _ fu,d U,d fu,d U,d _ fu,d U,d 

Cq Vq -I- X2 Vq X3 fg 



fu.d U,d I fu,d U,d , fu,d U,d I fu,d U,( 

Co ^0 +^1 ^1 +^2 ^0 +^3 ^3 



(54) 



where fu = {u, u}, fd = {d, I}, couple to the appropriate Higgs VEVs vf''^, respectively. 



Finally, the eigenstates corresponding to the eigenvalues Eq. (|51|) or Eq. (|53|) may be expressed 
as linear combinations of the fields defined in Table IIIH which are localized at the intersections 
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between D6 branes. One finds that the eigenfunctions assume the forms 





f -ax 1 


, IA4) oc 








' 1 




^ X I 







where X = ( I ; (55) 



where a — )• 1 in the special case of Eq. (j53p . Most significantly, the ratios of the third- to fourth- 
generation coefficients defining these linear combinations are both given by the same ratio of the 
Yukawa couplings on the second two-torus, yi/y2- Thus, so long as the Yukawa couplings yi and 
y2 are the same for both the up- and down-type fermions, no mixing occurs between the third and 
fourth generations at the trilinear order. This condition can be easily satisfied if the D-branes are 
split only on the first and third two-tori.^ 



B. Analysis 

A full comparison of model predictions to experimental data would require the full RGE evolu- 
tion of the observed masses and CKM elements from the weak scale to the unification scale Mgut- 
However, the masses of the fourth-generation fermions, and their mixings with the other three, are 



of course unknown 



3- 



As discussed above, the combination of direct observation bounds, elec- 



troweak (EW) precision tests, and perturbative unitarity constrain ts p lace strong constraints on 
the possible masses, which is just as true in the SM as in the MSSM [47']. However, four-generation 
models, both super symmetric and not, are plagued by the presence of several large Yukawa cou- 
plings that exhibit Landau poles at the TeV scale and above. While attempts have been made to 
stabilize the numerical evolution of the RGEs in four-generation models up to Mgut by including 
new matter fields {e.g., ^8]), questions of the robustness of such models remain. For the following 
analysis, we assume that the effective running of the Yukawa couplings between the GUT scale and 
the EW scale may be stabilized by some mechanism, and we note that vectorlike matter similar 
to that introduced in 43] is present in the model, as seen in Table HIl For simplicity, we further 
(crudely) assume that the running of the Yukawa couplings between the GUT and EW scales varies 
slowly, so that weak-scale values of observables may be used. 

Let us recall that the Pati-Salam gauge symmetry has been broken to the SM by brane splitting. 
For the moment, assume that this splitting occurs only on the first two-torus. In fact, to forbid 
this mixing just for quarks, ci = C2 on the second torus is sufficient. Then one can write the total 



Likewise, the two massless-state eigenvectors depend only upon the third-torus parameter combination zi/z2. 
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brane shifts on the first torus for each type of fermion, given by Eq. (j36p with r = 1 as 

V _ 2ea2 + — fcl u _ ^ggl + gfe ~ ^cl /^.^n 

I _ 2ea2 + Cfe — €c2 d _ 2€ai + efo — 

' ~ 4 ' ' ~ 4 

where quarks begin on stack al and leptons on stack a2, while the neutrinos and up- type quarks 
both end on stack cl, and the charged leptons and down-type quarks both end on stack c2. One 
then observes that 

e"-e'^ = ^(-eci + ee2) = e"-e'. (57) 

Taking, for example, e'^ = 0, one must then satisfy the constraint 

e'=' = e" + e^ (58) 

when choosing parameters to fit the mass ratios of third- and fourth-generation fermions. Note 
that although we are treating the brane-shifts as free parameters, these are open-string moduli that 
should be stabilized in a completely realistic compactification. This requirement can be satisfied 
if the D-branes wrap rigid cycles, which are available on several different backgrounds 



49|-l5lll. 



In addition, the Kahler moduli on each torus may 



3e stabilized by including supergravity fluxes 



3 



and/or by gaugino condensation in a hidden sector 

From the relations given in the previous subsection, we make certain general observations. First, 
Eq. (|54p shows that a very light third-generation neutrino Vr can readily be obtained if 

(x[; + X2>o"-xK-^3^3 «0. (59) 



where we take vq = V2 so that the generic eigenvalues of the Yukawa matrices are given by Eq. (I53p . 
This condition can easily be satisfied if, for example, the Higgs VEVs take the values 

= 2CI{xl + xl), vl = Clx\, vl = Cjxl (60) 

where C is constant. Using these values, Eq. ()53p shows that niy^ vanishes while niyi may be large. 
With these VEV choices, the neutrino Yukawa matrix takes the form 

{ YZ Yz\ 

3^" ~ 2C , (61) 

\YZ YZ j 

and essentially degenerates to rank 1, so that only v' receives mass at leading order. Note that 
this rank-1 condition is independent of the specific values assigned to the Kahler modulus and shift 
parameter . 
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0.40-1 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 1: Ratio of the t to t' quark masses as a function of the shift parameter e = e", where the Higgs VEVs 
are chosen to satisfy Eq. (|M)) . and k ~ ^/3n. 

Since the up-type quark masses depend upon the same Higgs VEVs as the neutrinos, the 
question now becomes whether it is possible to obtain masses for the t and t' quarks in the range 

nif' , , 

— 2-3 . (62 

At first glance, simultaneously satisfying these criteria may appear difficult. However, the Yukawa 
couplings x" are in general different than xj' since the position on the first two-torus of the al 
stack of branes upon which the quarks begin may be different than the position the al stack from 
which the leptons begin, so that ea\ 7^ ea2i and — = \[ia\ — ^ai) ■ 

We now show that one can obtain results consistent with Eq. ()62p for generic nonzero values 
of K = — iK^^) and e". For example, let us choose e*^ = and k = -v/Stt (an arbitrary value used 
purely for illustration), while the Higgs VEVs have the values given by Eq. (j60p . The resulting 
ratio mt/rriti given by Eq. (j60p as a function of the shift parameter e = e" is plotted in Fig. [TJ This 
plot shows that mass ratios satisfying Eq. (f62ll can be obtained for 

e"* 0.19-0.31, and 0.69-0.81 . (63) 



Similarly, the down-type quark masses and charged leptons depend upon the Higgs VEVs v^. 
However, in addition the shift parameters must satisfy Eq. (1580 . The Higgs VEVs may be chosen 
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FIG. 2: Mass ratio between the third- and fourth-generation down-type quark (and charged-lepton masses) 
as functions of Kahler modulus k — ^/Stt on the first torus and shift parameter e= e'^. 

freely, apart from overall normalization; for simplicity, let us choose 

v^, = 2D/ix"',+x%), vi = D/x'\, vi = D/x'l (64) 

where D is a constant and x'^ are the same functions used in calculating the eigenvalues in 
the neutrino sector, but with a different value of the input parameter k, which we label k' . Note 
that the Yukawa couplings for the down-type quarks and charged leptons must be calculated using 
the same Kahler modulus k = -v/Svr as used for the up-type quarks and neutrinos. However, the 
down- type Higgs VEVs are chosen by the ansatz Eq. (164p . with k' < k. Obviously, many other 
suitable choices are possible. 

In Fig. [2] we plot the ratio of third- to fourth-generation masses for the down-type quarks and 
charged leptons using n = \/37r and k' = \/2vr/2 (again, purely an illustrative choice), as a function 
of brane shift parameter £ = 6^. As one sees, the resulting function exhibits several maxima, with 
minima that are equally spaced about the maxima. The appearance of these alternating maxima 
and minima can, of course, be traced to the fact that the mass eigenvalues are essentially linear 
superpositions of the 'd functions exhibiting well-defined symmetry properties in e. 

In fitting the third- to fourth-generation mass ratios for the down-type quarks and charged 
leptons, we require 

~ 80-130 , ^ « 60-700 , (65) 

mft rrir 
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These mass ratios are clearly much larger than rriti/mt. Furthermore, one must impose Eq. (j58p 
on the shift parameters. Fortunately, as noted previously, the mass ratio between third- and 
fourth-generation down-type quarks and neutrinos plotted in Fig. [2] exhibit minima near e = 0.375, 
e = 0.125, which are separated by Ae ~ 0.25. By choosing values near e'^ ^ 0.375 and e' f« 0.625, 
satisfying all constraints is possible. Indeed, if one chooses = 0.3737 and = 0.1237, for 
example, one finds 



■my 
nib 



97.7 and 



267.3 . 



(66) 



Choosing nearly symmetric shift values naturally provides nearly equal Yukawa couplings for the 
b and r, so that 6-r unification may be naturally achieved. Since e^ = e'^ — = 0.25, one obtains a 
value of near a maximum of Fig [TJ One then obtains the phenomenologically acceptable value 

2.62 . (67) 

Furthermore, noting that 



mt 



XqVq + X^V^ + X2VQ + Xgfg 

where xj' and xf differ only by the use of the shift parameters e'^ = and e" = 0.25, respectively, 
one finds 



mt' 



0.924 . 



(69) 



Using the tabulated masses mt^b,T 



46j and these ratios, one obtains mf = 452.5 GeV, nib' = 
409.5 GeV mf-nib' = 43.0 GeV, and m^/ =474.9 GeV, m^/ = 418.1 GeV ^ m^/-m^/ = 56.8 GeV, 
meaning that all the mass constraints listed in the Introduction are satisfied. 

In addition, no mixing occurs between the third- and fourth-generation fermions since, as dis- 
cussed above, the D-branes have only been split on the first two-torus. We also note that one can 
split the D-branes on the third two-torus as well, without inducing mixing between the third and 
fourth generations. Doing so does not affect the mass ratios m^/ni^ for each type / of fermion. 
However, splitting on the third two-torus does affect the relative mass hierarchy between the up- 
type quarks and the neutrinos, as well as the between the down-type quarks and charged leptons. 
The suppression of fourth-generation quark mixin g im plies long quark lifetimes and the possibility 



of hadronic bound states 



52l | , as discussed in [53|, |5j] 



As seen in this analysis, one can obtain mass splittings within each generation in a fairly natural 
way, since the mass ratios involve superpositions of i? functions that have convenient symmetry 
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properties. As a result, the mass ratios between the third- and fourth-generation fermions exhibit 
maxima and minima as a function of shift parameter e. Note that generic points reflect fully rank-2 
Yukawa matrices, while the minima correspond [according to Eq. (|54p ] to points where the Yukawa 
matrices degenerate to rank 1. This tendency toward Yukawa matrices of lower rank was also 
observed in the rank-4 model of [15]. The mass ratio mt/mt' is constrained to be much larger than 
the ratio between mh/mi,', mr/mr', and especially my^/rriyi. However, one sees that these ratios 
can be obtained if the shift parameter is chosen to lie near a maximum of the mass ratio curve, 
while the shift parameters e'^, e', and correspond to minima. For the present model, the brane 
shift parameters e", e*^, e\ and are continuous parameters constrained by Eq. (j57p . 

IV. CONCLUSION 

We have presented a globally consistent four-family MSSM constructed from intersecting D6 
branes wrapping cycles on a T^/(Z2 x Z2) orientifold. In contrast to other similar constructions, 
the Higgs fields appear in the chiral sector rather than as vectorlike matter that occurs between 
nonintersecting D branes. We find the gauge couplings to be unified at the string scale and the 
hidden-sector gauge groups to become confining in the range lO^'^-lO^^ GeV. Thus, the model is 
a four-generation MSSM at low energies. 

The Yukawa matrices of the model are rank 2, so that only the third- and fourth-generation 
fermions receive masses at trilinear order. The first two generations, in principle, receive masses via 
higher-order corrections, and so should be naturally lighter. This result contrasts sharply with that 
of our earlier rank-4 model jisl , in which only small regions of the VEV parameter space accommo- 
date the known fermion mass and mixing spectrum. In addition, the third- and fourth-generation 
fermions have nondegenerate masses but do not mix at leading order. Finally, a numerical analysis 
of the mass ratios between the third and fourth generations shows that simple, elegant, and fairly 
natural solutions satisfying all constraints are relatively easy to find. In particular, we have shown 
that one can obtain three light neutrinos and a heavy fourth, while simultaneously obtaining heavy 
t and t' quarks in a completely natural way. In addition, one may easily obtain masses for the 
h' quark and r' lepton that are heavy in comparison to the b quark and r lepton, respectively. 
Moreover, 6-r unification may be naturally realized. We conclude that our analysis demonstrates 
that the putative fourth-generation fermion masses required to satisfy the constraints placed upon 
them by direct experimental observations, precision electroweak measurements, and perturbative 
unitarity constraints can emerge in a completely natural fashion. The very real possibility that a 
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fourth generation of chiral fermions exists deserves serious consideration, one that will be settled 
conclusively at the LHC. 
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